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Abstract. We present preliminary results of 2D radiative modellin g of He 1 
lines in solar prominences, using a new numerical code developed by us (|Leger et all 
120071 ) . It treats self-consistently the radiation transfer and the non-LTE statisti- 
cal equilibrium of H and, in a second stage, the one of He using a detailed atomic 
model. Preliminary comparisons with new visible plus near-infrared observations 
made at high spectral resolution with THeMIS are very satisfactory. 



1. Introduction 

Multi-dimensional radiative transfer is of general interest in various domains of 
astrophysics. Hereafter our primary interest concerns the radiative modelling of 
isolated and illuminated structures in which non-LTE plasma conditions prevail: 
solar prominences. In particular, spectral lines of He I are commonly used for the 
diagnosis of m agnetic fields in solar promi nences. While most inversion codes, 



diagnosis or magneti c fields m solar promi nences. Wnile most inversion codes, 
except hazel ( Asensio Ramos et al.ll2008l ). assume that spectral lines such as 



D3 and 1083 n m for instance are optically th in, it is not always confirmed obser- 
vationally (e.g. lLopez Ariste Sz Casmill2002l ). We also have a good idea about 



conditions under which g eometry effect s may have an impact on the formation 
of moderately thick lines (|Paletoulll997l ). In order to improve our ability to diag- 
nose prominences using such spectrosc opic data, we hav e thus developed a new 
2D numerical radiative transfer code ( Leger et al. 20071 ). Hereafter we present 



the numerical strategies we adopted and comment on the performances of this 
code. Finally, preliminary results of the 2D radiative modelling of He 1 lines and 
a first comparison with new observations are shown. 



2. Context and state-of-the-art 



Solar promi nences are dense and coo l structures hanging in a hot and low den- 
sity corona ( Tandberg-Hanssenlll995l ). They are often the source of CMEs, and 
the magnetic field very likely plays a major role in the triggering of those insta- 
bilities leading to the ons et of such large solar plasma ejections (see for instance 
Patsourakos Vial 200 ll . and references therein). 

He 1 lines such as D3 and 1083 nm are among the best tools for the study of 
magnetic fields in solar prominences. Full-St okes observations of these spectral 



lines, made at THeMIS (jPaletou et al.1 120011) have also led to a recent revision 



of magnetic field inversion tools ( Lopez Ariste Casini 2002). 
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Figure 1. Maximum relative change R c on the populations from an iteration 
to another against the number of iterations for the MG, SOR, Gauss-Seidel 
(GS) and ALI numerical schemes. MG provides a considerable gain in number 
of iterations and computing time. 



The most recent radiative models ( Labrosse k. Gouttebrozd l200ll . |2004| ) 
assume ID static slabs and no atomic fine structure for the He model-atom, 
leading to non-realistic gaussian synthetic triplet lines profiles. However this 
fine structure is conspicuous in our observations, and measurement of the ratio 
between the red and blue peaks of D% and 1083 nm are often in contradic- 
tion with the commonly assu med hypothesis of optically thin spectral lines (see 
Lopez Ariste &i Casinil 20021 ). 



It is therefore important to use the best numerical radiative transfer meth- 
ods, in 2D geometry at least, and a detailed Hei atomic model in order to 
improve our ability: i) to compute realistic synthetic spectral line profiles and 
then ii) to improve further available magnetic field inversion tools. 



3. Numerical strategies 



We adopted the following numerical strategies. The formal solution of the ra- 
di ative transfer equation is solved in 2D using the short characteristics method 
of [A ucr & Paletoul (|!994l ). To improve convergence, we used a Gauss-Seidel it- 
erative scheme with Successive Over- Relaxation GS/SOR for mu ltilevel atoms 
( Trujillo Bueno Fabiani Bendichol Il995l ; iPaletou Lege r 2007). Finally we 
included this 2D-GS/SOR iterative scheme into a nested multigrid method (see 
Fabiani Bendicho et al.l[l997l : iLeger et al.ll2007f ). 

Then, in order to synthesize He line profiles, we proceed as follows. We solve 
the statistical equilibrium and the radiative transfer equations for the hydrogen 
atom. Then the statistical equilibrium equations for the helium atom are solved 
consistently with the ionization equilibrium already computed for hydrogen. For 
the first step, we considered a plasma composed of neutral and ionized hydrogen, 
and neutral helium with an abundance by number fixed at 0.1. 
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Table 1. Computing time and number of iterations (between parenthesis) 
for four different numerical methods and two spatial grids. 



Grid size 


MALI 


GSM 


SOR 


MG 


163x163 


184mn29s (156) 


120mn27s (123) 


78mn41s (84) 


31mn46s (29) 


243x243 


698mn38s (259) 


417mn24s (201) 


269mn21s (130) 


65mn31s (41) 



As an illustrative example, we adopted a temperature T = 8000 K, a gas 
pressure p g = 0.05dyncm~ 2 and a microturbulent velocity £ = 5kms~ 1 . The 
geometrical width of the slab was D y = 5 000 km and its height D z = 30 000 km, 
sampled with a 243x243 points spatial grid. We considered just 5 bound levels 
and the continuum for hydrogen. 

The respective rates of convergence for the ALI, Gauss-Seidel (GS), SOR 
and multigrid (MG) multilevel iterative processes in 2D geometry are displayed 
in Fig. [TJ where it is plotted the maximum relative change R c on the populations 
from an iteration to another, against the number of iterations. In Tab. [H we 
also provide the corresponding computing time and number of iterations for two 
different spatial grids. 

This shows that the MG scheme is definitely superior, both in iteration 
numbers and computing time, against all other numerical schemes. This is 
especially true when the grid is more and more refined. And we ought to remind 
again here that grid refinement is very impo rtant for the sake of precision, as it 
was demonstrated by IChevallier et al. (|2003l ) . 



4. 2D synthetic line profiles vs. observations. 

We anticipate that more realistic models should also consid er the multi-th read 



or spatial fine structure, as observed in prominences (see iHeinzell [2007, and 
references therein). 

Hereafter, we consider a single thread whose geometrical width and height 
are taken identical, D VtZ = 1 200 km, and sampled with a 123x123 grid. We also 
adopted five different temperatures T =10 000, 12 000, 15 000, 16 000 and 17000 
K, a gas pressure p g = l.Odyncm -2 and a microturbulent velocity £ = 5kms~ 1 . 

Concerning the atomic models, we considered 5 bound levels and the con- 
tinuum for hydrogen, but 17 bound levels (up to n=3) and the continuum for 
helium, including the atomic fine structure for levels 2 3 P, 3 3 P and 3 3 D. We took 
explicitly into account the hydrogen pumping for UV l ines and the c ontinuum 



of He, and the coronal illumination for UV lines of He ( Tobiska 199 ll ). 



Synthetic profiles and a sample of our observations made at THeMIS in 
June 2007 are plotted in Fig. [2j At first glance, this preliminary comparison 
with observations is very promising since we are able to recover very nicely the 
high-resolution features of the line profiles. 
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Figure 2. Normalized (a) D3 and (b) 1083 nm spectral lines observed at 
THeMIS in June 2007. Theoretical normalized line profiles for (c) D3 and (d) 
1083 nm for 5 different temperatures are displayed for comparison. 



5. Conclusions 

We have briefly presented preliminary results of the 2D radiative modelling of 
He I lines in solar prominences, using a new 2D numerical radiative transfer code 
developed by us. So far, we favour 2D multi-thread models in order to explain a 
number of properties of observed He 1 spectral lines. A more detailed comparison 
with our new observations made at THeMIS, simultaneously in the visible and 
the near-IR, has begun. A next step will be to evaluate the impact of such a 2D 
radiative modelling on magnetic field inversion tools. 
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